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Aims. The P1/P2 period ratio of transversal loop oscillations is currently used for the diagnostics of longitudinal structuring of 
coronal loops as its deviation from 2 is intrinsically connected to the density scale-height along coronal loops and/or the sub-resolution 
structure of the magnetic field. The same technique can be applied not only to coronal structures, but also to other oscillating magnetic 
structures. 

Methods. The oscillations in magnetic structures are described by differential equations whose coefficients depend on the longitudinal 
structure of the plasma. Using a variational principle written for the transversal component of the velocity vector, developed earlier by 
McEwan et al. (2008), we investigate how the different temperature of the environment compared to the temperature of the magnetic 
m,j , structure will influence the P1/P2 ratio for typical coronal and prominence conditions. The possible changes are translated into 

. quantities that are used in the process of remote plasma diagnostics in the solar atmosphere. 

' Results. Using a straightforward, yet comprehensive, procedure we show that under coronal conditions the effect of the temperature 

(— I ' difference between the plasma inside and outside the magnetic structure can change considerably the period ratio; in the case of 

^ ^, coronal loops the change in the period ratio can reach even 40%. We also show that once dispersive effects are taken into account, 

I . with oscillation periods shorter than the cut-off period (determined by the density and temperature difference) the domain where the 

O ' model can be applied is reduced. In the case of prominences embedded in the hot corona, the effect of the environment is negligeable 

given the high density and temperature difference between the chromospheric prominences and corona. Using a numerical approach, 
C/3 , we show that our analytical discussion produces a robust result. We also discuss what implications our model has on seismological 

■ (or diagnostics) techniques in the solar corona. 

Conclusions. Our analysis shows that the period ratio P1/P2 is sensitive to the temperature difference between the loop and its 
environment and this effect should always be taken into account when estimating the degree of density structuring with period ratio 
^ , method of coronal loops. 

0^ , Key words. Magnetohydrodynamics (MHD) — Sun: corona — Sun: magnetic fields — Sun: oscillations 
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1. Introduction where observations of wave characteristics (amplitude, wave- 

length, propagation speed, damping time/length) are corrobo- 

^ Sudden energy releases m the solar atmosphere ai-e known to ^^^^^ ^j^^ theoretical modelling (dispersion and evolutionary 

^ . generate large scale global waves propagatmg over long dis- equations, as well as MHD models) in order to derive quantities 

. tances (see, e.g. Moreton and Ramsey 1960, Uchida 1970, j^at cannot be direcdy or indirectly measured (magnetic field 

^ ; Thompson et al. 1999, Ballai et al. 2005). The energy stored magnitude and sub-resolution structuring, ti-ansport coefficients, 

k> in these waves can be released by traditional dissipative mech- ^^^^^^^ functions, thermal state of the plasma, stratification pa- 

^ _ anisms, but it could also be transferred to magnetic structures r^meters, etc.). Considerable advances have been achieved in 

H , which may come in contact with global waves. This scenario diagnosing the state of the field and plasma (see, e.g.Roberts 

. ^ . IS true not only for coronal structures, but applies to all mag- 1934^ Nakariakov et al. 1999, Nakariakov and Ofman 

netic entities in the solar atmosphere that can serve as waveg- 20OI, Ofman and Thompson 2002, Ruderman and Roberts 2002, 

uides (see, e.g. WiUs-Davey and Thompson 1999, Patsourakos ^^dries et al. 2005, 2009, Ballai et al. 2005, 2011, Gruszecki 

and Vourhdas 2009, Liu et al. 2010, etc). In the corona, BIT 20O6, 2007, 2008, Banerjee et al. 2007, Ofman 2007, 

waves generated by coronal mass ejections (CMEs) and/or flai-es 2009, McLaughlin and Ofman 2008, Verth et al. 2007, Ballai 

could interact with coronal loops, resulting in the generation of 2007, Ruderman et al. 2008, Van Doorsselaere et al. 2008, Verth 

kink modes, i.e. oscillations which exhibit periodic movement g^^^jyi 20O8, Verth et al. 2008, Morton and Erdeiyi 2009, 

about the loop's symmetry axis. Global waves can also inter- Ruderman and Erdeiyi 2009, Andries et al. 2009, Selwa and 

act with prominence fibrils as observed by, e.g. Ramsey and Qfman 2009, Selwa et al. 2010). It is highly likely that higher 

Smith (1966), and more recendy by Eto et al. (2002), Jing et al. resolution observations made possible recendy by space satel- 

(2003), Okamoto et al. (2004), Isobe and Tnpathi (2007), Pinter jj^es such as STEREO, Hinode, SDO (and future missions) will 

et al. (2008). Oscillations of magnetic structures were (and ai-e f^^^ther help understanding the complicated reality of the solar 

currendy) used as a basic ingredient in one of solar physics ^^^^^^ environment. Indeed since their launch, data provided by 
most dynamically expanding fields, namely coronal seismology. 
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these satellites are already shedding light on numerous aspects of 
coronal seismology, e.g. Verwichte et al. (2009) used STEREO 
data to determine the three-dimensional geometry of the loop, 
SDO/AIA data was used by Aschwanden and Schrijver (2011) 
to prove the coupling of the kink mode and cross-sectional os- 
cillations that could be explained as a consequence of the loop 
length variation in the vertical polarization mode. Finally, based 
on Hinode data, Ofman and Wang (2008) provided the first evi- 
dence for transverse waves in coronal multithreaded loops with 
cool plasma ejected from the chromosphere flowing along the 
threads. On the other hand the development of even the funda- 
mental mode is not always guaranteed, as was shown observa- 
tionally by, e.g. Aschnwanden et al. (2002) and later using MHD 
modelling by Selwa and Ofman (2010) and Selwa et al. (2011 
a,b). 

The dispersion relations for many simple (and some quite 
complicated) plasma waves under the assumptions of ideal mag- 
netohydrodynamics (MHD) are well known; they were derived 
long before accurate EUV observations were available (see, e.g. 
Edwin and Roberts 1983, Roberts et al. 1984) using simpli- 
fied models within the framework of ideal and linear MHD. 
Although the realistic interpretation of many observations is 
made difficult especially by the poor spatial resolution of present 
satellites not being quite sufficient, considerable amount of in- 
formation about the thermodynamical and dynamical state of the 
plasma, and the structure and magnitude of the coronal magnetic 
field, can still be obtained. 

The mathematical description of waves and oscillations in 
solar structures is, in general, given by equations whose coef- 
ficients vary in space and time. It has been recognised recently 
by, e.g. Andries et al. (2005) that the longitudinal stratification 
(i.e. along the longitudinal symmetry axis of the tube that coin- 
cides with the direction of the magnetic field) modifies the pe- 
riods of oscillations of coronal loops. Accordingly, in the case 
of kink waves, these authors showed that the ratio P\/P2 (where 
Pi refers to the period of the fundamental transversal oscilla- 
tion, while P2 stands for the period of the first harmonic of the 
same oscillation) can differ - sometimes considerably - from the 
value of 2, that would be recovered if the loops were homoge- 
neous. These authors also showed that the deviation of P1/P2 
from 2 is proportional to the degree of stratification (see also, 
e.g. McEwan et al. 2006, Van Doorsselaere et al. 2007, Ballai 
et al. 2011). Later, studies by, e.g. Verth et al. (2007), showed 
that it is not only density stratification that is able to modify the 
P1/P2 period ratio, but the variation of the loop's cross section 
area has also an effect on the period ratio. While the density 
stratification tends to decrease the period ratio, a modification 
of the cross section (i.e. when the magnetic field is expanding 
as we approach the apex) tends to increase the P\/P2 value. 
Observationally is much easier to detect the fundamental mode 
of kink oscillations (having the largest amplitude and the small- 
est damping rate), however high resolution observations made 
possible the evidence of even higher harmonics (see, e.g. De 
Moortel and Brady 2007, Van Doorsselaere et al. 2009). 

In this study we will restrict our attention only to the effect of 
density stratification, the effect of the magnetic field structuring 
is left for a later analysis. The period ratio P1/P2 is connected 
to the density scale-height that quantifies the variation of den- 
sity along the magnetic structure. All previous studies consid- 
ered that the density stratification (indirectly the scale-height) is 
identical inside and outside the magnetic structure. However, the 
scale-height is directly linked to the temperature (via the sound 
speed) and an equal scale-height would mean an equal tempera- 



ture, clearly not applicable for, e.g. coronal loops and/or promi- 
nence fibrils. 

The aim of this paper is to investigate, using a simple math- 
ematical approach, the effect of the environment on the period 
ratio P1/P2 and the consequences of the inclusion of a distinct 
environment on estimations of the degree of density stratifica- 
tion. The paper is structured in the following way: in Section 2 
we introduce the mathematical formalism and obtain analytical 
results for typical coronal and prominence conditions. Later, in 
Section 3 we use our findings to draw conclusions on the im- 
plications on coronal seismology and we present a method that 
could help diagnosing not only the degree of stratification in the 
loop, but also the temperature ratio of the plasma inside and out- 
side the loop. Finally, our results are summarised in the last sec- 
tion. 



2. The mathematical formulation of the problem and 
analytical results 

EUV observations made by the recent high resolution space 
satellites (SOHO, TRACE, STEREO, Hinode, SDO) showed 
that, after all, coronal loops are enhancements of plasma (trac- 
ing the magnetic field in virtue of the frozen-in theorem) and 
the density of a typical loop can be as much as 10 times larger 
than the density of the environment. The heating of these coro- 
nal structures - according to the accepted theories (see, e.g. 
Klimchuk 2006, Erdelyi and Ballai 2007 and references therein)- 
occurs at the footpoints, while thermal conduction, flows, waves, 
instabilities and turbulences will help the heat to propagate along 
the full length of the loop. It is also obvious (as see in X-ray by, 
e.g. Hinode/XRT) that the temperature of the loop exceeds the 
temperature of the environment. A typical length of a coronal 
loop is 20-200 Mm, which means that the density inside the loop 
(seen in EUV) can vary by an order of magnitude, leading to the 
necessity of studying the effect of density stratification on the 
oscillations of coronal loops. 

As pointed out for the first time by Dymova and Ruderman 
(2005, 2006), the propagation of kink waves in a straight tube in 
the thin tube approximation can be described by 



0, 



(1) 



where v,- denotes the radial (transversal) component of the veloc- 
ity vector and the quantity ck is the propagation speed of kink 
waves (often called the density weighted Alfven wave) defined 
as (see e.g. Edwin and Roberts 1983) 



ck = 



2 2 \I/2 

Pi + Pe , 



(2) 



where p, and denote the densities inside and outside the coro- 
nal loop, and vai and vac represent the propagation speeds of 
the internal and external Alfven waves, respectively and the dy- 
namics is treated in the cold plasma approximation. For identical 
magnetic field inside and outside the tube, the fact that p, > p^ 
means that vas > vai- If we lift the thin flux tube restriction, then 
Eq. ([TJ must be complemented by terms that would describe dis- 
persion. It is interesting to note that a similar equation was found 
recently by Murawski and Musielak (2010) describing Alfven 
waves. As we specified, our approach is using the cold plasma 
approximation in which the dynamics of kink oscillations in a 
coronal loop is described by Eq. ([TJ. For the sake of complete- 
ness we need to mention that the cold plasma approximations 
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is not always true, especially in hot coronal loops observed by 
the SUMER instrument (e.g. Wang et al. 2003). It is very likely 
that the dynamics of kink oscillation will be described by a sim- 
ilar equations as Eq. ^ but with an extra term resulting from 
considering the effect of pressure perturbation. Furthermore the 
consideration of higher values of plasma-y6 will affect the values 
of eigenfrequencies as was found by McLaughlin and Ofman 
(2004), and Ofman (2010, 201 la). 

Assuming that all temporal changes occur with the same fre- 
quency, CO, we can consider that the temporal dependence of vari- 
ables (including v, ) has the form exp(/a;f), which means that the 
PDE given by Eq. ([T]i transforms into 



Vr 



0. 



(3) 



In reality the kink speed does not depend only on the longitudi- 
nal coordinate, z, but on all 3 coordinates. It is known that the 
dependence on the transversal coordinate, r, leads to the phe- 
nomena requiring short transversal length scales (resonant ab- 
sorption, phase mixing, turbulence, wave leakage), used to ex- 
plain the rapid damping of kink oscillations (see, e.g. Ofman 
and Aschwanden 2002, Ruderman and Roberts 2002, Ruderman 
2008, etc). Equation (O implies that the eigenfunctions, v^, are 
driven by particular forms of cxiz), through the particular profile 
of the quantities that make up the kink speed (density, magnetic 
field). Inspired from the eigenvalue problem of Rayleigh-Ritz 
procedure, McEwan et al. (2008) used a variational principle that 
allows the calculation of eigenvalues, w, - a method that is em- 
ployed by our analysis. Let us multiply the above equation by vv 
and integrate from the apex to the footpoint of the loop as 



v,—dz + u>' -f 
Jo dz^ Jo cU 



-dz = 0. 



(4) 



Using integration by parts in the first integral (taking into ac- 
count that for the fundamental mode v, (L) - dvr{Q)/dz = and 
for the first harmonic Vr(Q) = v^(L) = 0), the above equation 
simplifies to 



Jo c\{z) Jo \dz] 



(5) 



which results into the equation derived earlier by McEwan et al. 
(2008) 



2 

where 



(6) 



-dz. 



In order to express the eigenvalue of such problem, we consider 
some trial functions for v,. that satisfy the boundary conditions 
imposed at the footpoints and the apex of the loop. Since we 
are interested only in the characteristics of fundamental mode of 
kink oscillations and its first harmonic, we will assume that v, (z) 
will be proportional to cos(7rz/2L) for the fundamental mode and 
sin(7rz/L) for the first harmonic. It is obvious that these choices 
for eigenfunctions correspond to the homogeneous plasma, how- 
ever - as we show in the Appendix - the corrections to the eigen- 
function due to density stratification are rather small. 

The problem of how the kink speed depends on the longitudi- 
nal coordinate, z, is a rather delicate problem and only simplified 



cases can be solved analytically. For simplicity, let us consider 
that the magnetic field inside and outside of the coronal loop are 
identical and homogeneous, while the density varies exponen- 
tially according to 

/ 2L nz 
Pi(z) = Pi(0) exp — sm — 
\nHi 2L 



Pe(z) = Pe(0) exp I sin ^ 



(7) 



where p,(0) and pdO) are the densities inside and outside the 
loop at z = 0, i.e. at the the loop apex and //, and are the 
density scale-heights inside and outside the loop. Obviously the 
choice of density reflects a simplified description of the coronal 
loop model where plasma is isothermal and other further effects 
are neglected, however, this density profile allows us to obtain 
analytical results. A realistic description would require taking 
into account that the plasma is not isothermal (inside and outside 
the loop), the loop is curved and the density can depend on other 
coordinates, as well. This form of density dependence on the z 
coordinate was earlier used by, e.g. Verth et al. 2007, McEwan 
et al. 2008, Morton and Erdelyi 2009, Morton and Ruderman 
20 11, Morton et al. 20 11. With our chosen density profiles, the 
kink speed given by Eq. (|2]i becomes 



2Bl 



^APi(z) + Peiz)] 



2vi.(0) 



sin f ) + r' exp (;;|^ sin f ) ' 



(8) 



where Bq is the magnitude of the magnetic field, VAiiO) is the 
Alfven speed at the apex of the loop, and ^ is the density ra- 
tio, i.e. p,(0)/pe(0). Since the density outside the coronal loop is 
smaller than inside, we will consider that ^ > 1 . The quantities 
Hj and H,, are the density scale-heights and they are proportional 
to the temperature of the plasma. Here we denoted = H^IHj. 
Since the temperature of the loop is higher than its environment, 
we will take x ^ 1, so that the value x - ^ corresponds to 
an identical density variation with height inside and outside the 
loop and identical temperatures, x °° resulting in a constant 
density in the environment of the loop, while the limit x ~^ ^ 
represents a case when the plasma inside the loop is homoge- 
neous. 

Using the particular form of v, for the fundamental mode and 
its first harmonic, we obtain that in the case of the fundamental 
mode 



8L 



^2 ~ 



8yyi,(0) 



X /i(2y) + Li(2y) + rV 



(9) 



where we introduced the dimensionless variable y - LlnHi, 
lyix) is the modified Bessel function of order v, Ly{x) is the mod- 
ified Struve function of order v, and the index / stands for the 
fundamental mode. For the first harmonic we obtain that 

*! = ^. 

' 2L 
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Fig. 1. The variation of the P1/P2 period ratio with the temper- 
ature parameter, and the ratio L/tt//, for the case of a typical 
coronal loop (the density ratio, ^, is 2). 



Fig. 2. The relative variation of the P1/P2 period ratio with the 
temperature parameter, x, and the ratio L/H for the case of a 
typical coronal loop (the density ratio, ^, is taken to be 2). 



= TTTTT^ [Id2y) + Li(2y)] - 3 [/„(2y)+] 



+Lo(2y) + + 

TT 



■'W7)-(?))^S 



(10) 



where the superscript 1 in the expressions of *P[ and stands 
for the first harmonic. Now using Eq. Q for both modes we 
obtain that 



Pi 
P2 



^2' 



(11) 



Inspecting the above relations we can see that the period ratio 
P1/P2 does not depend on Alfven speed or loop length (they 
cancel out when calculating Eq.fTTli. For coronal conditions we 
plot the period ratio given by Eq. ( fTTT i for ^ = 2 with the vari- 
able y varying between and 10, although the larger values of y 
are rather unrealistic since y - 10 would correspond to a scale- 
height of 30 times shorter than the loop length (the scale-height 
corresponding to a typical temperature of 1 MK is 47 Mm). 
Another variable in our problem is the ratio of scale-heights (i.e. 
temperatures), sox will be varied in the interval to 1. 

The dependence of the P\/P2 period ratio on x and the ratio 
L/nHj for coronal conditions (^ = 2) is shown in Fig. 1 with 
the case discussed earlier by, e.g. Andries et al. (2005) corre- 
sponding to the value = 1. In addition to the ratio L/nHj our 
model prescribes a possible diagnostic of the temperature dif- 
ference between the loop and its environment. The importance 
of changes when the different temperature of the environment 
is taken into account can be shown in a relative percentage plot 
shown in Figure 2. The relative change was calculated as the 
percentage change of the results of our investigation compared 
to the case when x - 1 . As we can see, the changes in the 
domain corresponding to values of x close to 1 are not signif- 
icant. However, as the temperature of the environment becomes 
lower than the temperature inside the loop, this difference shows 
changes of the order of 10-20 % for values of x of up to 0.5, 



while for the cases with x near zero, the difference can be even 
40 % (for;t' - 0-2 and L/nHi = 2). Since the relative change is 
negative, it means that for the same value of Pi /f 2 calculated as- 
suming the same temperature the ratio, L/nHi is overestimated. 
A change of 25% in the period ratio occurring at approximative 
values of L/nHj - 0.8 and;^' - 0-65 would mean that for envi- 
ronment temperature that is 35 % less than the loop temperature, 
the scale-height is underestimated by about 25%. It is important 
to note that the density ratio, ^, does have an important effect 
of the variation of period ratio. An increase of ^ to the value of 
10 would result in relative percentage change reduction and the 
maximum value of the change is attaining its maximum value at 
33 % (for X = 0.2 and L/nHj = 1). 

The same analysis was repeated for prominence structures. 
These structures are known to be of chromospheric origin and 
show rather long stability. Prominence fibrils are surrounded by 
much hotter and less denser corona. For these structure we sup- 
pose that the density of the prominence two orders of magnitude 
times higher, i.e. we take ^ - 100. The typical temperature of 
prominences varies between 5x10^ and 10^ K, while the tem- 
perature of the surrounding corona can be even two orders of 
magnitude higher That is why, the value ofx is chosen to change 
in the interval 50-150. As we can see in Fig. 3, the changes of the 
period ratio P1/P2 for prominences does not show large varia- 
tion with X and an analysis of the relative change (compared 
to the case corresponding to x - 50) would reveal that these 
changes are of the order of 0. 1 %. 

Strictly speaking the form of equilibrium densities given by 
Eq. ^ is obtained after imposing an equilibrium of forces along 
the vertical direction (considering that the loop is vertical) when 
the forces created by pressure gradients are balanced by gravi- 
tational forces. Moreover, the density scale-heights given before 
are connected to the gravitational acceleration. In an isothermal 
plasma the density scale-height is given as 

4 

yg 

where cs is the sound speed and y is the adiabatic index. In 
this case, the effect of the environment is described by a similar 
equation as given by Eq. ([T]l, however, the governing equation is 
supplemented by an extra term that describes dispersive effects. 
Recently, Ballai et al. (2008) studied the nature of forced kink 
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„(50, 1 SO) 




Fig. 3. The same as in Fig. 1 but we plot the variation of P1/P2 
for prominences where ^ = 100 and;^^ varies between 50 and 
150. 



case of a homogeneous loop (understood in a local sense) where 
the propagation of kink oscillations is described by a Klein- 
Gordon equation, we can easily find that the period ratio is given 
by 



Pi 

Pi 



1 



-1/2 



(17) 



where the second term in the square bracket gives the deviation 
of f 1/P2 from 2. 

It is well known (see, e.g. Rae & Roberts 1982, Ballai et al. 
2006) that only those waves will be able to propagate in such 
magnetic structure whose frequencies are larger than the cut-off" 
frequency given by Eq. ( fTSI l. This condition will impose an up- 
per boundary on the applicability region of the variables L/ttH, 
and X- The solution of Eq. (fTSl l is a transversal kink oscillation 
propagating with speed ck followed by a wake that is oscillating 
with the frequency ojc- Due to the height dependence of den- 
sities, the cut-off" frequency will also depend on z- In terms of 
the variables used in our analysis, the cut-off frequency can be 
written as 



oscillations in a coronal loop and they obtained that the dynam- 
ics of the oscillations is given by an inhomogeneous equation of 
the form 



d^Vr Pi-pedVr 2 / x ^^^r «r 
or Pi + Pe OZ OZ-^ 



(12) 



where !F represents the external driver. This equation can be cast 
into a Klein-Gordon equation after introducing a new function, 
such that v,.(z, f) = Q{z, Oe''*'^-', where the value of A{z) is chosen 
in such a way that all first derivatives with respect to z vanish. 
Introducing this ansatz into Eq. ( fT2b we obtain 



5^ 



Pi - Pe 



Pi + Pe 
g(Pi-Pe) 



2cl{z)T{z) 



dQ 
dz 



Pi + Pe 



T{z)-ci{z)nz)' -ci{z)T{zY 



Q- 



&2 



Te 



(13) 



The condition that the coefficient of dQjdz is zero reduces to 

g(Pi - Pe) 



T{z) ^A'z + A 



^clizXpi + Pe) ' 



(14) 



With these restrictions, our governing equation transforms into 



&2 



Here tOc is the cut-off frequency of kink oscillations and is given 
by 

Mf:)- 



g^(pi-pe? _ gc\{z) d 

44(z)(p/+P.)2 2 dz 



Pi -Pe 



CK(z)iPi +Pe) 



(16) 



It is obvious that the existence of this dispersive term is due to 
the different densities between the loop and its environment. The 
dispersion will aff'ect the value of the P1/P2 period ratio. In the 



In 



2(£i+^-i£2) 



(^£i-r'£2)' 



vl(0) 



-—cos —(£1 £2) 

L 2L X 



1/2 



(18) 



where we used the notations 



exp 



£2 



exp^-sm- 



In order to evaluate the eff'ect of the cut-off" on the period ratio, 
P\/P2, let us now return to Eq. (fT2T l and repeat the same calcula- 
tion as before. For simplicity we neglect the inhomogeneous part 
on the RHS of Eq. ( fT2l i and assume that perturbations oscillate 
with the the same real frequency, oj. Using the same variational 
method, for the fundamental mode we obtain 



ip/ 1 VI//' 

2 ^1 ^1 



(19) 



given by 



Here and are already defined by Eq. ^ and T^' is simply 



8v%{0)L [ 



{2y) + Li(2y)- 



2(^-1) 



(15) 



(20) 



For the first harmonic we obtain that the frequency is given by 

(21) 



Vf/l + \J/1' 



VJ/1 

^2 



where and ^e specified by Eq. ( fTOl l and *Pj' is defined as 



2 

^ g 



1 



1 - -TT77^ ^o(2y) + Lo(2y) - -{h{2y) + U(2y)) 



'»'7)-'(f)-fH!)-(f 



(22) 
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P 1 /P2 for D 



= 0.05 & X — (50,150) 




0.8 l^-t^+X>"--"^ 1.2 




Fig. 4. The variation of P1/P2 period ratio for coronal loops 
when the effect of the cut-off period is taken into account . 



Fig. 6. The same as Figure 4, but now the period ratio is plotted 
for prominence conditions. 
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Fig. 5. The domain where the periods of fundamental mode and 
its first harmonic are smaller than the kink cut-off period, given 
by Eq. (fTSI l for coronal conditions. The domain of permitted val- 
ues is shown by label "I", while region "11" corresponds to the 
unphysical results. 



The influence of the kink cut-off period on the period ratio P1/P2 
is shown in Figs. 4 and 5 and its effect becomes obvious when 
Figs. 1 and 4 are compared. First of all, we need to note that the 
effect of the dispersive term is more accentuated for increasing 
values of L/nHj. Next, imposing the condition that the periods, 
we are investigating, are smaller than the cut-off period means 
that the domain of interest is restricted, as shown in Fig. 4 (in 
fact, only Pi is required to be smaller than the cut-off period, as 
f 2 is always smaller than Pi). The region where the above condi- 
tion is not satisfied was flagged by zero and the drop in the P1/P2 
to zero represents the boundary of the region where this imposed 
condition is satisfied. Looking from above, the domain of appli- 
cability is shown in Fig. 5, with the domain labelled by index 
"I" indicating the region where the periods are smaller than the 
cut-off period, while the region "11" corresponds to the set of 
(L/7THi,x) values for which no physical solution is found. When 
the dispersive term is not taken into account, the P1/P2 period 
ratio is independent on the length of the loop and the Alfven 
speed (basically, the independent parameter is the density ratio. 



f). However, once the dispersive term is considered in Eq. (fTZb . 
the P1/P2 ratio will depend on the length and the Alfven speed 
measured at the apex of the loop. For the numerical example 
shown here, we have chosen a length of 150 Mm and an Alfven 
speed at the apex of the loop (z = 0) of 1000 km s"'. Since the 
length of the loop is given, varying L/nHi would mean a change 
in Hi. In our numerical analysis we stopped at L/irHi = 5 that 
corresponds to a scale-height of 9.5 Mm. Assuming a loop in 
hydrostatic equilibrium, the value of H = 9.5 Mm leads to a 
temperature of 0.2 MK. As the value of L/jtHi decreases, the 
scale-height increases. It is easy to verify that for a fixed value 
of L, the P1/P2 period ratio is proportional to VAi(Q), and an in- 
crease/decrease of 200 km s~' would result in a change of only 4 
% towards the large L/nHj part of our investigated domain. If we 
fix the value of the Alfven speed at the loop apex, the variation 
of P1/P2 is inversely proportional to L, but again, for a change 
of 50 Mm in L drives changes of the order of 2% for large values 
of L/nHi. 

For solar prominences we repeated the calculations but now 
assuming that ^ - 100, L - I Mm and v^iiQ) - 120 km s"' (see 
Fig. 6). Under these conditions it is obvious that possible solu- 
tions are found for L/nHi < 2.4. The condition that a hydrostatic 
equilibrium is reached inside the prominence means that the 
smallest scale-height we use is 1.32 x 10^ m, which corresponds 
to a minimum temperature of 2800 K (nearly a quarter of the 
typical prominence temperature). As we approach smaller values 
for L/nHi, the temperatures increase, so that when L/nHi = 0.7, 
the temperature is approximately lO"* K, a typical prominence 
temperature. 

3. Implications for magneto-seismology 

The immediate implication of our calculations is that the P1/P2 
period ratio has no one-to-one correspondence with the internal 
stratification of the magnetic structure, but depends also on the 
temperature ratio between the interior and exterior of the mag- 
netic structure, i.e. an observed period ratio allows the diagnos- 
tics of the temperature ratio, too. Our analysis shows that the 
effect of temperature difference is more pronounced for those 
cases where the temperature inside the waveguide is larger than 
outside (e.g. coronal loops) and, in general, negligible in promi- 
nence cases. For coronal loops it is also evident that noticeable 
effects of the temperature difference on the P1/P2 ratio are en- 
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Fig. 7. An example on how the period ratio of the first three har- 
monics of a coronal loop kink oscillations can be used to diag- 
nose the density scale-height of the loop and the temperature dif- 
ference between the coronal loop and its environment. Here the 
P1/P2 dependence is shown by the dotted line while the solid 
line stands for the value of f i/f 3. 



countered for relative small values of L/ttHi (say, below 5) and 
temperature ratio that are smaller than 70%. 

Our analysis also opens a new way of diagnosing the multi- 
temperature loops and their environment. The relations derived 
in the present study show that the physical parameters entering 
the problem are the density ratio, temperature ratio and the ra- 
tio of loop length and scale-height. Out of these quantities, the 
density is the parameter that can be determined (although with 
errors) from emission, so we will suppose that the value of ^ is 
known. The diagnostics of the loop in the light of the new intro- 
duced parameter becomes possible once we specify an additional 
relation connecting the temperature ratio and the density scale- 
height measured against the length of the loop. As a possibility 
we investigate the case when for the same loop we can determine 
now only the period of the fundamental mode (Pi) and its first 
harmonic (P2) but also the period of the second harmonic, here 
denoted by P3. Now we can form a new ratio. Pi /P3, which can 
be determined in a similar way as above. Since the measurement 
of the three periods refer to the same loop we can estimate the 
value of X and L/nHj in a very easy way. In Fig. 7 we illustrate 
such a case. We suppose that for a loop of half-length of 150 Mm 
with an Alfven speed at the apex of the loop of 1000 km s ' we 
measure the period ratio of 1.72 for P1/P2 and 2.67 for P1/P3. 
By specifying the value of the period measurement means that 
in a dependence similar to the one shown in Fig. 1 we obtain 
an arc which in a {LlnHi,x) coordinate system will look like 
the dotted line curve in Fig. 8. In a similar way, for the value of 
P1/P3 = 2.67 we would obtain another curve, here shown by 
the solid line (here the curves are the projections of the intersec- 
tion of similar surfaces as in Fig 1 with the horizontal surface 
corresponding to the specified period ratio). Since the two mea- 
surements correspond to the same loop, their intersection point 
will give us the exact value of L/tt//, and;^'. For the particular 
example used here we obtain LjTiHi - 0.57 wd x - 0.67. Our 
results show very little sensitivity with the density ratio, for ex- 



ample, if the ratio would be 10 then the intersection point would 
change to LjnHi - 0.54 anAx - 0.66. 



4. Conclusions 

In order to carry out coronal seismology it is imperative to know 
the relationship between the composition of a plasma structure 
and the oscillations supported by the coronal loop. High reso- 
lution observations make possible an accurate diagnostic of not 
only the magnetic field strength, but also the thermodynamical 
state of the plasma. 

The period ratio P1/P2 (and its deviation for the canonical 
value of 2) belonging to the period of transversal fundamental 
kink mode and its first harmonic is a perfect tool for diagnosing 
the longitudinal structure of magnetic structures. In our study 
we investigated the effect of the environment on the period ratio 
assuming that the density scale-heights (implicitly the tempera- 
ture) inside and outside magnetic structures are different. 

Using a simple variational method first applied in this con- 
text by McEwan et al. (2008), we derived for the first time an 
analytical expression that connects the value of the period of 
kink oscillations and parameters of the loop. We showed that 
in the case of coronal loops the effect of temperature difference 
between the loop interior and exterior can lead to changes of 
the order of 30-40% that could have significant implications on 
the diagnosis of longitudinal density structuring of the coronal 
loop. In the case of prominences, due to the very large density 
and temperature difference between the prominence and coronal 
plasma, the changes in P1/P2 due to the different temperature 
are very small. Once dispersive effects are taken into account 
(through a Klein-Gordon equation) the domain of applicability 
of P1/P2 seismology in the case of coronal loops becomes re- 
stricted and physically accepted solutions can not be found for 
any temperature ratio (here denoted by the parameter^). 

Since our model introduces a new variable in the process 
of plasma diagnostic a new relation is needed that connects the 
parameters of interest. To illustrate the possibilities hidden in our 
analysis we have chosen the case when the same loop shows the 
presence of the additional second harmonic. Superimposing the 
dependences of the P1/P2 and P1/P3 ratios with respect to the 
temperature ratio factor, x and LlnHi we could find the set of the 
values that satisfies a hypothetical measurement. 

Finally we need to emphasise that our approach supposes a 
certain degree of simplification, therefore our results do not pro- 
vide an absolute qualitative and quantitative conclusion. First we 
supposed that the loop is thin, and Eq. ([T) can be applied to de- 
scribe the dynamics of kink oscillations in coronal loops. It is 
obvious that this statement is not true for very short loops (the 
ratio of the loop ratio and its length is not very small) in which 
case, the governing equation has to be supplemented by an extra 
term. Secondly, our isothermal supposition of the loop and its en- 
vironment is also that needs refinement as observations (see, e.g. 
Winebarger et al. 2003, WaiTen et al. 2008, Berger et al. 2011, 
Mulu-Moore 201 1) show that the loops are not always in hydro- 
static equilibrium nor isothermal. Here we supposed the idealis- 
tic situation of a static background, however recent analysis by 
Ruderman (2011) showed that the the temporal dependence of 
density through flow and cooling can also influence the ratio of 
the two periods. 
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Appendix A: Corrections to the eigenfunctions due 
to the density stratification 

In the Appendix we estimate the corrections to the chosen eigen- 
functions due to the density stratification. Analytical progress 
can be made in the small ylx limit. Since is a value smaller 
than one, this condition would automatically mean that we work 
in the small y limit, provided x is not becoming too small. We 
are interested only in the characteristics of fundamental mode 
of kink oscillations and its first harmonic. Following equation 
Q with the boundary conditions vv(i) = dvriQ)/dz - and 
Vr(0) - VriL) = for the fundamental mode and first harmonic, 
we introduce a new variable so that 

^.£, with 

In the new notations the density inside the loop can be written as 



PiiO = Pi(0) exp 



h nC 
— sin — 

Hi 2 



(A.l) 



where h is the loop height above the solar atmosphere (a similar 
equation can be written for the external density). Next, we are 
working in the approximation /i/H, = e <s; 1, so equation ([3]l 
becomes 



1 + esin i —\ + - ( 1 + -sin i —\ 



v,(f) = 0,(A.2) 



where ^ - pj{0)/pe(Q) > 1 is the density ratio and;^^ = He I Hi < 
1 , with He and //, being the density scale heights inside and out- 
side the loop. Let us we write vv and O as 

v,(0 = vf\0 + ev;."(0 + Oie"), Q = Qo + eQi + 0{^).{K3) 



Substituting these expansions into Eq. ( IA.2b and collecting terms 
proportional to subsequent powers of e we obtain 



2„(0) 



0, 



(A.4) 



d^^ 



+ (1 + r')^yr^ + 2(1 + r')f^oniv<.°^+ 



(A.5) 



These equations must be solved separately for the fundamental 
mode and its first harmonic taking into account the boundary 
conditions 

- For the fundamental mode 

Vrii = 1) = 0, = 0) = 

d{ 

- For the first harmonic 

Vr(^ = 1) = y,(^ = 0) = 
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Let us first calculate the correction to the fundamental mode. It 
is easy to show that the solution of equation ( IA.4I ) taking into 
account the above boundary condition becomes 



and 



2vrTr 



(A.6) 



(A.7) 



It is important to note that the form of the above solution is ex- 
actly the same as the solution we employed for the eigenfunc- 
tion, V,-. In the next order of approximation we obtain equation 
(IA.5b which can be written as 



2 



2(1 + r')i^of^ivf ^ - fi^vr 1 1 + ^ " 



l2„(0) 



(!)■ 



(A.8) 



This boundary value problem will permit solutions only if the 
right-hand side satisfies the compatibility condition that can be 
obtained after multiplying the left-hand side by the expression 
of vf^ and integrating with respect to the variable ^ between 

d^^Q. 



and 1, or 
-1 



j| [2(1 + r')f^oOivf + Q^vf (1 + sin (^) 
After some straightforward calculus we can find that 



1 



1.1 



3(1+^-1)3/2 

As a result, equation ( IA.5I ) becomes 



(A.9) 



d^v'P 

1 v 

de 4 ' 



(1) 



(i+f) 



1 



■ cos 



(f) 



-^sin(^)cos(^) 



0. 



(A. 10) 



This differential equation will have the solution 
v^')(^) = C,sin(|) + C2Cos(f) 

+ ^) 2 cos (f ) + sin (f ) + "^sininO 



3n 



(1+^) 



(A. 11) 



Applying the boundary condition v|.''(^ = 1) = 0, we find the 
constant 



Ci =- 



(1^ 
3(1 +^)' 



In order to find the value of C2, we use the property of orthogo- 
nality, i.e. 



Jo 



vf\^^'> dC = 0, 



Correction to v_{r}(z) when e=0.1 



(%) 




1 ' 1 ' 1 ' 1 ' 1 ' 1 ' 1 ' r 

2. X 10'4. X 10''6. X 10''8. X lO^l.x 10* 1.4x10* 



Fig.A.l. Correction to the eigenfunction for the fundamental 
mode kink oscillation when e - 0.1. Here L - 1.5 x 10** m 
represents the loop length 

which result in 



C2 



9(1 H-^)7r 



As a result, the first order correction to Vr corresponding to the 
fundamental mode is 



(1+a 



1 . ln(\ 7 I7:^\ 

— sin — cos — 

3 I 2 / 97r I 2 / 



(1+^) 



l2cos(f) + .^sin(f) + ^sm(.0 



(A. 12) 



In Figure A.l we plot the correction to the eigenfunction for 
e - 0.1,^ - 0.9, and ^ - 10. Figure dA.lb shows that we can 
approximate Vr{z) by cos(7rz/2L) since the first order correction 
brings changes of about 1(%), i.e. insignificant. 

A. 1 . Corrections to the first iiarmonic 

The same analysis can be repeated for the first harmonic, taking 
into account the right boundary conditions. After a straightfor- 
ward calculation it is easy to show that 



(A. 13) 



1+^ 

16 ^+l/x 
I5n 1+^ 



176 . , 4 ^ 

sminO + 77 cos(;r^) 
2257r 15 



5;r /7r<r\ 37r /3<\ 

cos — H COS - 

32 I 2 / 32 \ 2 I 



Y cos(;rO - ^ sininO 



(A. 14) 
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T ' 1 ' 1 ' 1 ' 1 ' 1 ' 1 ' 1 

2. X lO'^ X 10''6. X lO'^S. X lO'^l.x 10* 1.4x10* 



Fig. A.2. The same as Figure A. 1 but here we represent the cor- 
rection to the eigenfunction for the first harmonic kink oscilla- 
tion. 



The correction to the eigenfunction corresponding to the first 
harmonic has been plotted in Figure A. 2 for the same values as 
before. It is obvious that the changes introduced by stratification 
in the value of the eigenfunction are of the order of 2(%), i.e. 
negligably small. The two figures show that the efi'ect of den- 
sity stratification becomes more important for higher harmonics. 
Given the very large values of x we used for prominences, the 
approximations used in this Appendix will always be vahd. Fo 
the graphical represention of corrections in Figs Al and A2 we 
used X - 1. If we lower this value to, e.g. 0.7 the corrections 
would still be small since the maximum relative change in the 
eigenfunction describing the fundamental mode would be I.l 
%, while for the first harmonic, this would increase to 3.5 %. 

The robustness of our analysis was checked using a full nu- 
merical investigation for arbitrary values of x and y. A typical 
dependence of the P1/P2 period ratio with respect to L/nHi for 
one value of x is shown in Fig. A. 3 where the solid line cor- 
responds to the analytical and the dotted line represent the nu- 
merical results, in both cases the density is inhomogeneous with 
respect to the coordinate z. The loop is set into motion using 
a gaussian-shaped source and we use a full reflective boundary 
conditions at the two footpoints of the loop. After the oscillations 
are formed, we use the FFT procedure to obtain the values of 
periods. Our analysis shows that the differences between the re- 
sults obtained using the variational method and a full numerical 
investigation are of the order of 7% but towards the large range 
of L/nHi. Restricting ourself to realistic values, i.e. L/nHi < 5, 
we see that the results obtained with the two methods coincide 
with great accuracy. 



